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1 Abstract 

^1 , The effect of the random magnetic field distribution on the phase diagrams and ground state mag- 
netizations of the Ising nanowire has been investigated with effective field theory with correlations. 
Gaussian distribution has been chosen as a random magnetic field distribution. The variation 
of the phase diagrams with that distribution parameters has been obtained and some interesting 
results have been found such as disappearance of the reentrant behavior and first order transitions 
which appear in the case of discrete distributions. Also for single and double Gaussian distributions 
, ground state magnetizations for different distribution parameters have been determined which 
can be regarded as separate partially ordered phases of the system. Keywords: Ising Nanowire; 
random magnetic field; Gaussian magnetic field distribution 

2 Introduction 

Recently there has been growing interest both theoretically and experimentally in the magnetic 
nanomaterials such as nanoparticles, nanorods, nanotubes and nanowires. Nowadays, fabrication of 
these nanomaterials is no longer difficult, since development of the experimental techniques permits 
us making materials with a few atoms. For instance, acicular magnetic nano elements were already 
fabricated [JJ [2J [3] and magnetization of the nanomaterial has been measured [4] . Nanoparticle 
systems have growing application areas, e.g. they can be used as sensors [5], permanent magnets [5J, 
beside some medical applications 7 . In particular, magnetic nanowires and nanotubes have many 
applications in nanotechnology [5J [5] ■ Nanowires can be used as an ultrahigh density magnetic 
recording media [TU1 HU H2] and they have potential applications in biotechnology [131 [H] , such as 
Ni nanowires can be used for bio seperation [151 [To] . 

In the nanometer scale, physical properties of these finite materials are different from those of 
their bulk counterparts. Some properties of these materials, which highly depend on the size and 
the dimensionality, can be used for fabrication of materials for various purposes. From this point 
of view, it is important to determine the properties of these materials theoretically. Most common 
used theoretical methods for determining the magnetic properties of these materials are mean field 
approximation (MFA), effective field theory (EFT) and Monte Carlo (MC) simulation as in bulk 
systems. For instance, nanoparticles investigated by EFT with correlations [T7], MFA and MC 
[18] . The phase diagrams and the magnetizations of the nanoparticle described by the transverse 
Ising model have been investigated by using MFA and EFT [Jj|J[5D]. Moreover, investigation of 
compensation temperature of the nanoparticle [21] and magnetic properties of the nanocube with 
MC [22] are among these studies. 

Another method, namely variational cumulant expansion (VCE) based on expanding the free 
energy in terms of the action up to m th order, has been applied to the magnetic superlattices [2"5] 
and ferromagnetic nanoparticles [241 2 51 . The first order expansion within this method gives the 
results of the MFA. 

Various nano structures can be modeled by core-shell models and these models can be solved 
also by MFA, EFT and MC such as FePt and Fe^O^ nanotubes [26,. The phase diagrams and 
magnetizations of the transverse Ising nanowire has been treated within MFA and EFT [27j [28] , 
the effect of the surface dilution on the magnetic properties of the cylindrical Ising nanowire and 
nanotube has been studied [29] [30], the magnetic propeties of nanotubes of different diameters, 
using armchair or zigzag edges has been investigated with MC [31 , initial susceptibility of the Ising 
nanotube and nanowire have been calculated within the EFT with correlations [32, 33 and the 
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compensation temperature which appears for negative core-shell coupling has been investigated 
by EFT for nanowire and nanotube [34] . There are also some works dealing with hysteresis 
characteristics of the cylindrical Ising nanowire [351 136] . Beside these, higher spin nanowire or 
nanotube systems have also been investigated, such as spin-1 nanotube [37] and nanowire [38], 
mixed spin - 3/2, 1 core shell structured nanoparticle [39], mixed spin - 1/2, 1 nanotube [40] systems. 

On the other hand, as far as we know, there have less attention paid on quenched randomness 
effects on these systems, except the site dilution. However, including quenched randomness or 
disorder effects in these systems may induce some beneficial results. For this purpose we inves- 
tigate the effects of the random magnetic field distributions on the phase diagrams of the Ising 
nanowire within this work. As stated in [29] the phase diagrams of the nanotube and nanowire 
are qualitatively similar, then investigation of the effect of the random magnetic field distribution 
on the nanowire will give hints about the effect of the same distribution on the phase diagrams of 
the nanotube. 

The Ising model in a quenched random field (RFIM) has been studied over three decades. 
The model which is actually based on the local fields acting on the lattice sites which are taken 
to be random according to a given probability distribution was introduced for the first time by 
Larkin [41] for superconductors and later generalized by Imry and Ma [42] . Beside the similarities 
between diluted antiferromagnets in a homogenous magnetic field and ferromagnetic systems in 
the presence of random fields (43J [44] , the importance of the random field distributions on these 
systems comes from the fact that, random distribution of the magnetic field drastically affects 
the phase diagrams of the system, and hence the magnetic properties. This situation has been 
investigated widely in the literature for the bulk Ising systems. For example, using a Gaussian 
probability distribution, Schneider and Pytte [45] have shown that phase diagrams of the model 
exhibit only second order phase transition properties. On the other hand, Aharony [46j and Mattis 
[47| have introduced bimodal and trimodal distributions, respectively, and they have reported the 
observation of tricritical behavior. With the same distributions and using EFT with correlations, 
Borges and Silva 48, 49, 50! showed that three dimensional lattices show tricritical behavior while 
two dimensional lattices do not exhibit this behavior. On the other hand, by using two site EFT 
instead of one site EFT, tricritical behavior can be observed on a square lattice [51) . Similarly, 
Sarmento and Kaneyoshi [52] investigated the phase diagrams of RFIM by means of EFT with 
correlations for a bimodal field distribution, and they concluded that reentrant behavior of second 
order is possible for a system with (q > 6). Recently, Fytas et al. [53] applied MC simulations on a 
simple cubic lattice. They found that the transition is continuous for a bimodal field distribution, 
while Hadjiagapiou [S3] observed reentrant behavior and confirmed the existence of a tricritical 
point for an asymmetric bimodal probability distribution within the MFA based on a Landau 
expansion. 

In a recent series of papers, phase transition properties of infinite dimensional RFIM with 
symmetric double [55] and triple |56j Gaussian random fields have also been studied by means of 
a replica method and a rich variety of phase diagrams have been presented. The situation has also 
been handled on 3D lattices with nearest-neighbor interactions by a variety of theoretical works 
such as EFT [57] [58], EFT with multi site spin correlations [59], MC simulations [60] ED 12] , pair 
approximation [63 , and the series expansion method[64). 

As seen in the short literature in bulk Ising systems, random field distributions keep up to date 
in the literature. Thus the aim of this work is to inspect the effects of random field distributions on 
the phase diagrams of the nanowire system as a nanostructure. The paper is organized as follows: 
In Sec. [3] we briefly present the model and formulation. The results and discussions are presented 
in Sec. 0] and finally Sec. [5] contains our conclusions. 

3 Model and Formulation 

We consider a nanowire which has geometry shown in Fig. [T] The Hamiltonian of the nanowire is 
given by 

H = -J\ ^2 S i S 3 ^ S m S n - J 3 22 S i S ™ ~ yHjSj - } H m S m (1) 

<i,j> <m,n> <i.m> i m 
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Figure 1: Schematic representation of a cylindrical nanowire (top view). The gray/black circles 
represent the surface/shell magnetic atoms, respectively. 



where s, is the z component of the spin at a lattice site i and it takes the values s,- = ±1 for 
the spin- 1/2 system. J\ and J 2 are the exchange interactions between spins which are located 
at the core and shell, respectively, and J 3 is the exchange interaction between the core and shell 
spins which are nearest neighbor to each other. Hi and H m are the external longitudinal magnetic 
fields at the lattice sites i and m respectively. Magnetic fields are distributed on the lattice sites 
according to a given probability distribution. The first three summations in Eq. (fT]) are over the 
nearest-neighbor pairs of spins, and the other summations are over all the lattice sites. 

This work -as a continuation of the earlier work [69j - deals with the following continuous 
magnetic field distribution, 

P (Hi) = P G (0, a) + ^ [G (H Q , a) + G (-H Q , a)} (2) 

where G (Hq, a) is the Gaussian distribution centered at Hq with a width a and it is given by 



1 \ V2 



G(H ,a) = - — - exp 
1 Zira^ ' 



(Hi - go) 
2a 2 



(3) 



Distribution given in Eq. ([2]) reduces to the system with zero magnetic field (pure system) for 
p = 1,(7 = 0. According to the distribution given in Eq. ([2]), p percentage of the lattice sites are 
subjected to a magnetic field chosen from the Gaussian distribution which has a width and Hq = 
as a center. Half of the remaining sites are under the influence of a field Hi which is randomly 
chosen from the magnetic field distribution G (Ho, a), whereas the distribution G (—Ho, a) used as 
distribution function on the remaining sites. 

Four different representative magnetizations (mi,i — 1,2,3,4) for the system can be given 
by usual EFT equations which are obtained by differential operator technique and decoupling 
approximation (DA) |65[ 166] . 

4 2 
mi — [Ai + miBi] [A 3 + m2B 3 } [A 3 + m 3 B 3 ] [A\ + m 4 i?i] 

m 2 = [A 3 + m 1 B 3 }[A 2 +m 2 B 2 } 2 [A 2 +m 3 B2] 2 

m 3 = [A 3 + m^g] 2 [A 2 + m 2 B 2 } 2 [A 2 + m 3 B 2 } 2 1 ' 

TO4 = [Ai + niiBi] 6 [Ai + m^Bi} 2 

Here mi, are the magnetizations of the two different representative sites in the core and m 2 , m 3 
are the magnetizations of the two different representative sites in the shell. The coefficients in the 
expanded form of Eq. Q are given by 

A\A\B™B n q = J dH t P (Hi) cosh fc ( J P V) cosh' ( J,V) sinh m ( J P V) sinh" ( J,V) / (H % , x) \ x=0 (5) 
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where V is the usual differential operator in the differential operator technique and the values of 
indices p, q can be p, q = 1, 2, 3. The function is defined by 



/ (Hi,x) = tanh (fix + 0Hi) . 



(6) 



as usual for the spin-1/2 system. In Eq. ([6]), /3 — l/(fcsT) where ks is Boltzmann constant and 
T is the temperature. The effect of the exponential differential operator to an arbitrary function 
F(x) is given by 

exp (aV)F 0) = F (x + a) (7) 

with any constant a. DA will give the results of the Zernike approximation |67j for this system. 
With the help of the Binomial expansion, Eq. (f4]) can be written in the form 



where 



mi 



m 2 



1 2 



m 3 2 77137714 



i=0 j=0fc=0i=0 
12 2 



m 3 = 



7774 



EEE^(M,fc)™i 

i=0j=0k=0 
2 2 2 

EEE^ 3 fc ) TO t 

i=0.j"=Ofc=0 
6 2 

EE^ 4 (i, ™>, 
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i=0 (=0 



(8) 



Kx (i,j,k,l) 
K 2 k) 
K 3 k) 
Ki (i, /) 




i-o-k A i-i B j+k B i 



A\~ l - l B\ +l 



(9) 



These coefficients can be calculated from the definitions given in Eq. ([5]) with using Eq. ([7]). 

For a given Hamiltonian and field distribution parameters, by determining the coefficients from 
Eq. ^ we can obtain a system of coupled non linear equations from Eq. ([5]l. and by solving this 
system we can get the magnetizations 777^,1 = 1,2,3,4. The magnetization of the core (m c ) and 
shell (m s ) of nanowire, as well as the total magnetization (ttit) can be calculated via 



(6777! + 7774) , 777 s 
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(6m 2 + 67773) , 777 T 
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(6777l + 6t77 2 + 67773 + 777 4 ) (10) 



Since in the vicinity of the critical point all magnetizations are close to zero, we can obtain 
another coupled equation system for determining this critical point by linearizing the equation 
system given in Eq. ([5]), i.e. 

A.m = 



where 



.4 



( #1(1,0,0,0) -1 #i(0,l,0,0) #i(0,0,l,0) #i(0,0,0,l) \ 

#2(1,0,0) #2(0, 1,0) -1 #2(0,0,1) 
#3(1,0,0) #3(0,1,0) # 3 (0,0,1)-1 
\ #4(1,0) # 4 (0,1)-1 / 

I 7771 \ 

m 2 
m 3 
\ 7774 J 



(11) 



(12) 



(13) 
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Critical temperature can be determined from det(A) =0. As discussed in [33j . the matrix A 
given in Eq. (|12[) is invariant under the transformation J3 — > — J3 then we can conclude that the 
system with ferromagnetic (J3 > 0) core-shell interaction has the same critical temperature as 
that of the system with anti- ferromagnetic ( J3 < 0) core-shell interaction (with the same | J3 1) for 
certain Hamiltonian and magnetic field distribution parameters. Although this discussion has been 
made for the system with zero magnetic field in [33 , this conclusion is also valid for this system, 
because of the symmetry of the magnetic field distribution. Equation det(A) = is invariant under 
the transformation J3 — » — J3 for the nanowire with magnetic field distribution given in Eq. ([5]). 

Beside this, there are some symmetry properties of the coefficients defined by Eq. (j9)). Although 
the numerical integration in Eq. (JS|) do not take too much time, using these symmetry properties 
will shorten the numerical calculation time. These symmetry properties can be found in Sec. [5] 



4 Results and Discussion 

In this section we discuss the effect of the continous random magnetic field distribution on the 
phase diagrams of the system. Since the phase diagrams are the same for the J3 > and J3 < 
with the same J3I, we focus ourselves on the case J3 > 0, i.e ferromagnetic core-shell interaction. 
We use the scaled interactions as 

J X = J, r n = ^, n = 2,3. (14) 
We start with a single Gaussian magnetic field distribution. 



4.1 Single Gaussian Distribution 

The form of single Gaussian distribution which is defined in Eq. @ is governed by only one 
parameter a , which is the width of the distribution. This distribution distributes negative and 
positive valued magnetic fields -which are chosen from the Gaussian distribution- to lattice sites 
so that sum of all lattice site's magnetic field is equal to zero. Although the total magnetic field 
is zero, randomly distributed negative and positive fields drag the system to the disordered phase. 
On the other hand, the interactions J n ,(n = 1,2,3) enforce the system to stay in the ordered 
phase. Another factor is the temperature which causes thermal agitations which induce disordered 
phase when energy supplied by the temperature to the system is high enough. Thus a competition 
takes place between these factors. 




Figure 2: The phase diagrams of the nanowire with single Gaussian random field distribution in 
the (ksTc/ J, a) plane for different r 2 ,r3 values. 



The phase diagrams in (fcsT c /J, a) plane can be seen in Fig. [2] for different r 2 and r 3 values. 
Firstly we can see from Fig. [2] that the general effect of increasing a is to decrease the critical 
temperature, as expected. There is not any evidence of reentrant behavior in the phase diagrams, 
as in ordinary Ising lattices 68 . The critical temperature value at a = and the critical value of 
the a -which makes the critical temperature zero- depend on r 2 ,r3, i.e. when r 2 or r% raises, these 
two critical values gets higher except one region. As in the same system with discrete distributions 
[69] when r 3 = 0.0 up to a certain value of the r 2 , the phase diagrams are the same in (/c_b T c /J, a) 
plane. In Fig. [2ja) all phase diagrams are the same for r 2 < 1.5. This fact can be explained by 
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taking into account the total spin-spin interaction strengths of the core and shell which do not 
interact with each other when 73 = 0.0. 

In the absence of the interaction between the core and shell, the critical temperature of the sys- 
tem is determined by the core, up to a certain r 2 > 1.0 value. Two distinct shell spins (which have 
magnetizations 7712,7713) have four nearest neighbors, and two distinct core spins (magnetizations 
of which are labeled as 7711,7714) have five and eight nearest neighbors, respectively. Thus at the 
value of r 2 — 1.0, core has higher critical temperature than the shell. As r 2 increases, the critical 
temperature of the shell gets higher and after a certain value of r 2 , the shell begins to determine 
the critical temperature of the system. Before this value, changing r 2 can not change the phase 
diagrams of the system. The dependence of the critical temperature values of the core and shell 
can be seen in variation of their magnetization curves with temperature which are given in Fig. [3] 
for some selected values of r 2 with fixed r 3 = 0.0 and a = 3.0 values. 

Within this region, although increasing r 2 does not change the critical temperature of the system 
for a fixed value of a, it causes change in the magnetization behavior with the temperature as seen 
in Fig. [SJc). As seen in Fig. [3f b) , as r 2 increases for a fixed a, the ground state magnetization 
and the critical temperature of the shell layer increase. When r 2 reaches the value that makes the 
critical temperature of the shell equal to the critical temperature of the core, increasing r 2 after 
this value manifests itself in rising critical temperature of the system. 

We can see from Fig. [2] that, increasing r 2 and 73 values make the ferromagnetic region wider 
in (fcsT c /J, a) plane (except for the region explained above), since increasing r 2 and 73 values 
raises the absolute value of the lattice energy coming from the spin-spin interaction, which must 
be overcomed by both thermal agitations and random magnetic fields for the phase transition from 
an ordered phase to a disordered one. 




Figure 3: Variation of the magnetization with temperature for some selected r 2 values of the 
nanowire with single Gaussian random field distribution. The fixed parameter values are r% — 

0. 0, cj — 3.0. 

As we done in our earlier work |69) for discrete distributions, now we want to investigate the 
variation of ground state magnetization behaviors with a. Then our question is; how the ground 
state magnetizations change as er is varied for different r 2 , 73 values? We know from the conclusions 
of our earlier work , bimodal and trimodal discrete random magnetic field distributions produces 
a number of ground state magnetizations which are different from zero (disordered phase) and one 
(ordered phase). These intermediate states may be regarded as partially ordered phases. In Fig. 
IH we plot the variation of the ground state magnetization values (mg, mf , rrij. core, shell and 
total ground state magnetization, respectively) as a function of the width of the single Gaussian 
distribution (a) for some selected values of r 2 with 73 = 1.0. We can see from the Fig. 2] that, 
this distribution creates higher number of these partially ordered phases, in comparison with the 
discrete distributions. Although increasing a decreases the critical temperature of the system 
gradually as seen in Fig. ©, EFT based on DA approximation gives that, increasing a decreases 
the ground state magnetization as a general trend, but it creates some plateaus i.e. there are 
some regions where the ground state magnetizations does not change with varying a values. With 
increasing temperature, the effect of the increasing a on the magnetization becomes continuous 

1. e. increasing a decreases the magnetization continuously at higher temperatures as seen again in 
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Fig. H 




Figure 4: Variation of the ground state magnetizations with a for some selected values of for 
the system with single Gaussian magnetic field distribution. Each magnetization calculated at a 
temperature fc^T/ J = 0.001 can be regarded as ground state of the system. At the same time, 
the magnetization values calculated at the temperature ksT/J — 0.1 are plotted as dashed lines 
with the same T2 values. Fixed parameter value is = 1.0. 



4.2 Double Gaussian Distribution 

Let us investigate the phase diagrams of the nanowire with double Gaussian distribution. This 
distribution is given by Eq. ^ with p = 0. In Fig. [SJ we plot the phase diagrams in the 
(ksT c / J, Ho/ J) plane for some selected r^-, r3 and a values. First, we see from the Fig. [SJthat, for 
fixed r2, t~3 and a values, when the distance between the centers of the Gaussians (2Hq/ J) increases, 
the critical temperature of the system gets smaller. Beside this, for fixed r2,r% values, increasing 
a makes the ferromagnetic region in the (fcsT c / J, Hq/ J) plane narrower, since the randomness 
increases due to increasing a. Another effect of the rising a in this distribution is that, increasing 
<7 destroys the reentrant behavior which exist in the system for given r2, values for the bimodal 
distribution (i.e. double Gaussian distribution with a — 0.0) as seen in Figs. [5ja),(d),(e),(f). Not 
only reentrant behavior disappears when a rises, but also first order transitions and tricritical 
points. At the same time as seen in Figs. Efb) and (c), rising a also destroys other first order 
transitions which originate from an ordered phase to a disordered phase and which may not come 
with the reentrant behavior. Then we can say that for the nanowire, while bimodal distribution can 
induce a first order transitions -and also tricrtical points- at higher Hq/ J values, double Gaussian 
distrbitubion can destroy these first order transitions when a is large enough. 

The ground state magnetizations appear in more complicated forms for this distribution than 
those corresponding to the bimodal distribution|69j . This can be seen from Fig. [6] which represents 
the ground state magnetizations corresponding to the phase diagrams given in Fig. [Sic) . As in the 
single Gaussian distribution, large number of partially ordered phases appear in the system which 
disappear when the temperature rises. 



5 Conclusion 

In this work, the effect of the continuous random magnetic field distributions on the phase diagrams 
and ground state magnetizations of the Ising nanowire has been investigated. 

There are two distribution parameters (Hq / J and a) of the random magnetic fields which can 
also be regarded as degree of the randomness which is imposed by the distribution of the magnetic 
field on the nanowire. Rising Hq/J or a can be considered as higher degree of randomness. This 
shows itself as a change in the phase diagrams which are plotted in (fcsT c / J —a) or (fcsT c / J — Hq / J) 
planes. We can say that, this change is two fold when we compare the results of the continuous 
distribution with the discrete distribution. First one is obvious; the ferromagnetic region in that 
planes collapses with rising degree of randomness. The second one is, rising degree of randomness 
imposed by rising a on the system, destroys all first order transition lines and also tricritical points 
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Figure 5: The phase diagrams of the nanowire with double Gaussian random field distribution in 
the {ksT c / J 1 Hq/ J) plane for different r2,r^^a values 



at which the first order and second order transition lines meet. Rising degree of randomness which 
is caused by rising er, also destroys reentrant behavior. While the system with magnetic field 
distribution given in Eq. ([2]) with a — 0.0 can show reentrant behavior, after a certain a -which 
depends on and ra- this behavior disappears. This means that, the system with double Gaussian 
distribution can not pass the border of the two phases with thermal agitations from a disordered 
phase to an ordered phase after a certain a value. 

Rising degree of randomness also shows itself in the ground state magnetizations. As one can 
see in |69j . bimodal magnetic field distribution can induce a few plateaus in the (rnL — KgT/J) 
plane. These plateaus where the ground state magnetization values do not change with rising 
randomness can be regarded as partially ordered phases. We can see from Figs. 0] and |5] that, 
continuous random magnetic field distribution can induce higher number of these partially ordered 
phases. Also it can be seen in Fig. [4] and |6] that, it is impossible to relate the magnetization values 
or widths of the constant ground state magnetization region of the partially ordered phases, with 
the system parameters. 

When the temperature varies a little, these partially ordered phases disappear. This means 
that there have to be some first order transition lines ending with isolated critical points in Figs. 
[2] and [H However in order to clarify whether this partially ordered phases are artifact or not of 
EFT based on DA, it is necessary to investigate the problem with more advanced methods. 

We hope that the results obtained in this work may be beneficial form both theoretical and 
experimental point of view. 
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Figure 6: Variation of the ground state magnetizations with a for some selected values of Ho/ J 
for the system with double Gaussian magnetic field distribution. Each magnetization calculated at 
a temperature ksT/J = 0.001 can be regarded as ground state of the system. At the same time, 
the magnetization values calculated at the temperature ksT/J = 0.1 are plotted as dashed lines 
with the same Hq/J values. Fixed parameter values are r 2 = 1.0 and r$ = 1.0. 



6 Appendix A: Symmetry properties of the coefficients 

As seen in the definitions of the the coefficients given by (J9j> , interchanging between certain indices 
do not change the part of that coefficient which is defined by (0. For obtaining the simple 
symmetry properties originated from this fact, let us start with the first coefficient by defining 



K[ (i,j,k,l) 



i -i 



B\ +l Bi 



(15) 



From the definition given in Eq. ([5]) we can see that Eq. (|15j) is invariant under the transformations 
i — > I, I — > i and j — > k,k —tj, thus we can use this property as 



Kx (i,j,k,l) 
K x (l,j,k,i) 
Kx (i,k,j,l) 
Kx (l,k,j,i) 




K[ (i,j,k,l) 
K[ (i,j,k,l) 
K[ (ij,k,l) 
K[ (ij,k,l) 



(16) 



i.e. 



Kx (l,j,k,i) 
Kx (i,k,j,l) 
Kx {l,k,j,i) 



Kx (i,j,k,l) 

Kx (i,j,k,l) 
l 



(17) 



Kx (i,j,k,l) 



Eq. ([T7|) reduces the total number of coefficients which has to be calculated as separately from 60 
to 45. In a similar way, for the other coefficients we can obtain similar symmetry properties as 



K 2 (i,k,j) = K 2 (i,j,k) 



(18) 



K 3 (i,k,j) = K 3 [i,j,k) 



Ki(l,i) 



K 4 (i,l) 



(19) 
(20) 
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These properties reduce the number of separate coefficients from 18 to 6, 27 to 9, 21 to 15, 
respectively. Thus for this system, total number of coefficients reduces from 126 to 75. 

On the other hand, another important symmetry property of the coefficients which is useful 
in the calculations of the integration given in Eq. ([5]) is about the magnetic field. Let us choose 
P (Hi) = 5 (Hi — Hq) as a magnetic field distribution and denote the coefficient defined by Eq. 
J5]) as Qkimn (Ho, Jp, Jq) = ApA l qB™Bq. By writing hypergeometric functions in terms of the 
exponentials then writing them with binomial distribution and using Eq. (J7J we get 



2 ) y^y^y^y^ Crsfo/ (h , a p j p + a q j q ) (21) 

where 



r=0s=0t=0v=0 



C rstv = (*]('](?](") (-ir + "-^ (22) 



r I \ s I V t J \ v 
and 

dp = 2r + It — m — k, a q = 2s + 2v — I — n (23) 

Now we can see from Eq. (|21[) that, for each term, C rs t v f (Ho,a p J p + a q J q ) has a corresponding 
term which is Ck- r i^ s>m -t n ^ v f(Ho,—a p J p — a q J q ) in the sum. It can be seen from Eq. (|22l) 
that, these two coefficient are related to each other by 

Ck—r,l—s,m—t,n—v = (~ 1) C rs tv (24) 

which means that for odd m+n, the terms which have / (Hq, a p J p + a q J q ) and / (Hq, —a p J p — a q J q ) 
are opposite signed coefficients. If to + n is even then the terms which have / (Ho,a p J p + a q J q ) 
and / (H , —a p J p — a q J q ) will be the same signed coefficients in the sum given in Eq. (|2"Tj) . Thus 
we can conclude that, 

n ijT T t \ / ©fcimn (H , J p> Jg) , m+n even , , 



Qkimn (H , J p , J q ) , m+n odd 

Since the functions defined in Eq. (j6]) have the property / (Hi,x) — —f(—Hi,—x), then we 
can write Eq. (|2"5)) as 

n { u t t \ / -®khnn(H ,J p ,Jq), m+n even , . 

0WmB Jp ' J?) = 1 e Hm „ (H 0) J p , J 9 ) , m+n odd ^ 

Since the distribution given in Eq. ([2]) is symmetric about Hq = 0, then in integration in Eq. 
(J5j> for any term Qkimn (Hq, J p , J q ) has a corresponding Qkimn (—Hq, J p , J q ) term. This means 
that in the case where random field distribution is given by Eq. ([2]), 

Qkimn (Jp, Jq) = J dH P (Hq) Qkimn (Hq, J p , Jq) I _^ ^ m+n odd 

(27) 



is valid. Eq. ([27|) also reduces the number of the coefficients which has to be calculated, approxi- 
mately half of the whole set of the coefficients. 
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